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Abstract 

The  aim  of  this  investigation  is  to  construct  an  adaptive  observer  and  an  adaptive  compensator 
for  a  class  of  infinite-dimensional  plants  having  a  known  exogenous  input  and  a  structured  pertur¬ 
bation  with  an  unknown  constant  parameter,  such  as  the  case  of  static  output  feedback  with  an 
unknown  gain.  The  adaptive  observer  uses  the  nominal  dynamics  of  the  unperturbed  plant  and  an 
adaptation  law  based  on  the  Lyapunov  redesign  method.  We  obtain  conditions  on  the  system  to  en¬ 
sure  uniform  boundedness  of  the  estimator  dynamics  and  the  parameter  estimates,  and  convergence 
of  the  estimator  error.  For  the  case  of  a  known  periodic  exogenous  input  we  design  an  adaptive 
compensator  which  forces  the  system  to  converge  to  a  unique  periodic  solution.  We  illustrate  our 
approach  with  a  delay  example  and  a  diffusion  example  for  which  we  obtain  convincing  numerical 
results. 


1  Introduction 

In  this  paper,  we  construct  adaptive  observers  for  the  following  class  of  infinite-dimensional 
systems  on  a  given  Hilbert  space  X  with  inner  product  and  norm  {•,  •)  and  |  •  \x ,  respectively: 

=  Apx(t)  +  Bu(t )  +  f(t)\  x(0)  =  xq  E  X  (1.1) 

y(t)  =  Cx{t ),  (1.2) 

where 

Ap  =  A0  +  BTC  (1.3) 

and  Ao  is  a  generator  of  an  exponentially  stable  Co-semigroup  T  it) ,  i.e.  there  exist  constants 
M,  n  >  0  such  that 

\\T(t)  ||  <  Me_/rf.  (1.4) 

The  signals  u(t)  and  y(t)  are  the  vector- valued  inputs  and  outputs,  respectively.  f(t )  is 
an  X- valued  known  exogenous  input  and  B  E  £(C  m.  X).  C  E  C(X.  4’  m).  Consequently, 
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Ap  generates  a  strongly  continuous  Co-semigroup  for  any  T  and  the  system  is  well-posed. 
A  key  assumption  is  that  the  original  system  (Aq,B,C)  satisfy  a  positive-real  condition, 
where  by  positive-real  we  mean  the  following. 

Definition  1.1  Suppose  that  (?(•)  :  (Fq  C{dJm),  where  (Fq"  =  {s  :  Res  >  0}.  If 

(*)  G(s )  =  G(a);  (1.5) 

(ii)  G(s)  is  holomorphic  on  (Fq  ;  (1.6) 

(in)  G(s )  +  G(s)*  >  0  for  all  s  =  joj,  lu  EM,  (1.7) 

then  G  is  positive- real. 

In  fact,  we  require  that  G(s  —  y)  —  C  ((s  —  n)I  —  Ao)  1  B  be  positive-real  for  some  y  >  0. 

The  motivation  is  that  the  perturbation  BTC  arises  as  the  unknown  constant  gain  in 
the  feedback  law  u  =  Yy.  Thus,  the  actual  formulation  is 

x(t)  —  A0x(t)  +  BYy(t)  +  Bu(t)  +  /(t);  a?(0)  =  x0-  (1.8) 

Some  preliminary  results  on  adaptive  observers  of  this  type  were  given  in  Demetriou  and 
Ito  [10],  but  they  only  applied  to  exponentially  stable,  dissipative  systems  (Ao  +  Aq  <  -2  III) 
under  a  co-location  assumption:  B  =  C* .  In  this  case,  G(s  —  y)  is  automatically  positive- 
real.  In  Curtain,  Demetriou  and  Ito  [6],  the  co- location  assumption  was  circumvented  by 
appealing  to  a  (very  special  version)  of  the  positive-real  lemma.  Unfortunately,  this  version 
only  applies  to  a  system  for  which  Aq  is  self-adjoint,  exponentially  stable,  and  the  inputs  and 
outputs  are  scalar-valued.  In  this  paper,  we  introduce  a  third  type  of  positive-real  lemma 
based  on  recent  results  in  Curtain  [5] .  Although  there  are  many  results  on  the  positive-real 
lemma  in  the  infinite-dimensional  literature  (see  [21], [22], [23], [24],  and  [25], [26], [27]),  most 
are  in  terms  of  a  certain  Riccati  equation.  For  our  application,  we  need  the  following  singular 
version,  for  which  no  corresponding  Riccati  equation  exists.  We  require  the  existence  of 
operators  Q  €  C(X)  and  L  :  X  — »  C  m  such  that  for  x  €  D(A) 

A*Qx  +  QAx  =  -L*Lx  (1.9) 

B*Qx  =  Cx.  (1.10) 

The  earlier  results  in  the  literature  are  too  restrictive  for  our  application.  Balakrishnan  [1] 
assumes  that  A  is  Riesz-spectral  and  that  the  scalar  inputs  and  outputs  are  very  smooth; 
both  Curtain  [3,  4]  and  Pandolfi  [17]  require  exact  controllability,  which  is  never  satisfied 
by  our  class  of  systems.  In  the  recent  results  by  Curtain  [5]  and  Pandolfi  [18],  the  lat¬ 
ter  assumption  is  removed  and  [5]  provides  the  type  of  positive-real  lemma  suited  to  our 
application. 

In  section  2,  we  summarize  and  compare  three  distinct  versions  of  the  positive-real 
lemma  from  [10],  [6]  and  [5].  In  [10]  and  [6]  the  L-operator  maps  from  X  to  X,  in  contrast 
to  the  situation  in  the  usual  positive-real  lemma,  where  L  maps  from  A'  to  (D  m.  We  discuss 
three  examples  which  satisfy  at  least  one  of  the  three  versions  of  the  positive-real  lemma. 

In  section  3,  we  propose  the  following  adaptive  observer  for  the  structurally  perturbed 
system  (1.9)  with  the  observation  (1.2)  under  the  assumption  that  G(s  —  y)  =  C((s  —  y)I  — 
Ao)~1B  is  positive-real  for  some  y  >  0. 

x(t)  =  A0x(t)  +  Bu(t)  +  BY  (t)y(t)  +  /(t);  x(0)  =  x0 

f  (t)  =  GCe(t)yT(t );  f  (0)  =  f0. 
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(1.11) 

(1.12) 


0  <  G  =  Gt  is  an  adaptation  matrix  gain  and  we  only  require  that  /  be  locally  Bochner 
integrable.  The  main  result  is  that  if  u,  y  E  1/00(0,00;®  m),  then  T(t)  and  x(t)  and  the 

jU  1  X 

observation  error  e(t)  =  x(t)  —  x(t)  are  bounded  in  norm  for  t  >  0  and  \\e^lQ^  (x{t)  - 
x(t))\\x  — >  0  as  t  -»  oo,  where  Q  is  the  operator  in  the  positive-real  lemma  (1.9)  with  A  = 
Ao+pl.  If  Q  is  invertible  or  if  y  E  LooDL2(0,  oo;  (D  m),  we  obtain  \\eP(x(t)  —  x(t))||x  — >  0  as 
t  ->  oo.  Under  a  persistently  exciting  type  condition  we  can  achieve  parameter  convergence: 
T(t)  — >  T  as  t  — >  oo. 

In  section  4,  we  propose  an  adaptive  compensator  design  using  the  separation  scheme 

u(t)  =  u2{t)  -  f  {t)y(t)  (1.13) 

and  using  an  LQR  design  on  the  resulting  adaptive  observer 

x(t)  =  A0x(t)  +  Bu2{t)  +  f(t).  (1-14) 

Assuming  a  known  periodic  exogenous  input  /,  we  obtain  an  optimal  controller  u2(t)  of  the 
form 

u2(t)  =  -R~1B*  (Px(t)+r(t)),  (1.15) 

where  P  and  r(t )  depend  on  the  LQR  design  parameters  and  on  f{i).  The  corresponding 
closed  loop  trajectory  x(t)  converges  exponentially  to  a  periodic  signal  p(t).  The  overall 
adaptive  compensator  (1.11)  -  (1-15)  makes  the  closed  loop  trajectory  of  (1.8)  converge  to 
a  unique  periodic  signal  p(t)  in  the  following  sense 

(x(t)  —  p(t)) \\x  — >  0  as  t  — >  oo 

for  some  (3  >  0.  If  Q  is  invertible,  x(i)  converges  exponentially  fast  to  pit). 

To  illustrate  the  above  results  we  present  some  numerical  results  on  our  three  examples 
in  section  5. 

2  Positive  Real  Lemmas 

The  adaptive  observer  scheme  is  only  applicable  to  positive-real  systems  and  the  key  is  a 
positive  real  lemma.  Just  as  in  finite  dimensions  (1.5)-(1.6),  in  infinite  dimensions  it  is 
possible  to  have  different  versions  corresponding  to  spectral  factors  of  different  dimensions. 
We  have  found  three  useful  versions.  The  first  version  is  particularly  useful  for  dissipative 
systems  with  co-located  actuators  and  sensors,  and  it  was  utilized  in  Demetriou  and  Ito 
[10].  These  systems  are  always  positive- real,  and  the  following  lemma  is  trivial. 

Lemma  2.1  Suppose  that  A  is  the  infinitesimal  generator  of  a  contraction  semigroup  on 
X  and  B  E  C((Um:  X) .  Then  Q  =  I  is  a  solution  to  the  constrained  Lyapunov  equation  for 
x  E  D(A) 

(Ax,  Qx)  +  ( Qx ,  Ax)  <  0 

B*Q  =  B*. 

In  the  adaptive  observer  application,  one  also  needs  to  suppose  that  A  generates  an  expo¬ 
nentially  stable  Co-semigroup.  An  example  satisfying  Lemma  2.1  is  the  following. 
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Example  2. A  Consider  the  diffusion  equation 


dz  d2z 
dt  dx2 


+  b(x)u(t ); 


y(t)  =  /  b(x)z(x,t)da 

Jo 


z(0,t)  =  0  =  z{l,t) 


where  b  E  1/2(0, 1)  =  X. 
We  let 

D(A)  =  { 


and  define 


h  E  £2(0, 1)  :  h,  —  are  absolutely  continuous, 
„  dx 

d2h 

2  e  £2(0, 1)  and  h( 0)  =  0  =  h(  1) 


d2h 


Ah  =  for  h  E  D(A). 


Then  A  has  compact  resolvent,  eigenvalues  \n  =  —n2 7r2,  n  E  IN  and  eigenvectors  <pn  = 
y/2  sin(n7nr),  n  E  K  .  which  form  an  orthonormal  basis  for  Lotf).  1).  A  is  exponentially 
stable,  self-adjoint  and  for  x  E  D(A) 


(x,  Ax)  <  —\\x 


and  it  generates  a  contraction  semigroup. 

Finally,  note  that  y(t)  =  ( b,z(-,t ))  =  z(-,t)  and  B  =  C*. 

To  treat  systems  for  which  the  actuators  and  sensors  were  not  co-located,  the  following 
version  was  proven  in  Curtain,  Demetriou  and  Ito  [6]. 


Lemma  2.2  Suppose  that  A  is  self-adjoint,  has  compact  resolvent,  its  eigenvalues  {A n,n  E 
N  }  are  simple  and  its  eigenvectors  {<fin,n  E  JN  }  form  an  orthonormal  basis  for  X .  Suppose 
that  b,c  E  X  satisfy 

(c,  4>n)  (b,  4>n)  >  0,  n  E  IN 


sup 

neZV 


(ci  <fn) 

{b,  <f>n ) 


<  OO. 


Then  there  exist  0  <  Q  =  Q*  E  C(X),  L  E  C(D(A),  X)  and  y>0  such  that  for  x,y  E  D(A) 


{{A  +  yI)x,Qy)  +  (Qx,(A  +  yl)y)  =  -(Lx,Ly) 

(x,c)  =  { x,Qb ). 

Proof  Show  by  direct  substitution  that  the  following  operators  satisfy  the  constrained 
Lyapunov  equation 


T  _  f  —  2 {h  +  A„)(c,  4>n)  ^ 

h\  Xh)  ) 


{ X ,  <f,n)4>n ■ 


(2.2) 


□ 
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This  lemma  applies  to  Example  2. A,  where  we  can  take  n  =  tt2  —  e  for  any  e  >  0. 

The  following  example  from  Curtain,  Demetriou  and  Ito  [6]  does  not  satisfy  the  condi¬ 
tions  of  Lemma  2.1,  but  Lemma  2.2  does  apply. 

Example  2.B  Consider  the  diffusion  equation 

dz  d  2z  dz 

m  =  w  -  aai  +  i(x)v(t>  + 1 


z(0,  t)  =  0,  z(l,t)  =  0; 


with  output  given  by 


y(t)  =  [  e 
Jo 


where 


b(x)  = 


cz(x,t)  dx , 

1  on  [0,  |) 


0  elsewhere 

V. 

Take  X  =  1/2(0, 1)  to  be  the  Hilbert  space  with  the  weighted  inner  product 


(/.S')  =  [  e 

Jo 


Then,  defining 


D(A0)  =  { 


dh 


cf(x)g(x)  dx. 


h  :  h.  —are  absolutely  continuous 

dx 

d2h 

and  €  X,  h( 0)  =  0  =  h(  1) 


and 


.  ,  d2h  dh  .  ,  „.  .  . 

A0h  =  -j^  ~  a—  for  h  G  D(A0), 


it  is  straightforward  to  show  that  Ao  is  self-adjoint  with  eigenvalues  Xn  =  —  —  n2 tt2  and 
normalized  eigenvectors  < pn(x )  =  y/2eax/2  sm(mrx),  n  G  IN  .  The  set  {4>m  n  G  K  }  forms 
an  orthonormal  basis  for  X.  Let 


Cn  ■ —  (c-i  (frn)  — 


4mr\/2  (l  —  e  2  (— l)n^j 
4n27r2  +  a2 


,  n  G  K 


bn  ■ —  {bi  <pn)  — 


4nvrV2  (l  +  e"t  cos(^)  -  sin(^f )) 


4n27r2  +  a2 


So  bncn  >  0  for  all  n  and  for  certain  constants  m  and  M 


m  <  sup 

n>  1 


<  M. 


So  the  assumptions  of  Lemma  2.2  are  satisfied,  Q  given  by  (2.1)  satisfies  the  constrained 
Lyapunov  equation  (2.2)  and  it  is  boundedly  invertible;  L  is  also  bounded. 

Note  that  in  both  Lemmas  2.1  and  2.2  the  L  term  will  be  unbounded  in  general,  even 
though  B  and  C  are  bounded,  and  that  L  maps  into  the  state-space  X.  This  is  in  contrast 
to  the  usual  finite-dimensional  version  for  which  L  maps  into  the  output  space  (D  .  The 
latter  version  is  much  harder  to  prove  for  infinite-dimensional  systems  and  earlier  versions 
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in  Balakrishnan  [1],  Curtain  [3,  4]  and  Pandolfi  [17]  assumed  very  strong  conditions  on  the 
system,  such  as  exact  controllability.  Here  we  extract  some  useful  results  from  Curtain  [5], 
who  assumes  only  mild  conditions  on  the  system  operators  (A,B,C). 

First  we  need  some  extra  notation: 

.  .  _  J  /  :  Cq"  — i >  Z  and  /  is  holomorphic  and 

°°(  )  =  j  ||/||oo  =  BUPwgR||/(jW)||jr<00 


H2(Z) 


f  :  (D  J  — t  Z  and  /  is  holomorphic  and  [ 

ll/lll  =  suPx>o/^o  ll/(®  +  <  oo  J 


L2{{~joo1  joo);  Z) 


f  :  (—joo,  joo)  — >  Z  and  /  is  measurable  and 
ll/lb  =  (f-o o  Wf(ju)Wxdu)  2  <  oo 


H00(Z)  is  a  Banach  space  under  the  sup  norm  and  H2{Z )  and  L2((— joo,  joo);  Z)  are  Hilbert 
spaces  under  their  ||  •  ||2  norms.  /  €  H2(Z)  has  a  limit  to  a  function  /  E  L2((0,  joo);  Z)  and 
/  is  isomorphic  to  /  with  ||/||2  =  ||/||2-  /  and  /  are  usually  identified  with  each  other  and 
with  this  identification  H2(Z )  is  a  subspace  of  L2  ((—joo,  joo);  Z).  We  denote  by  Ph2  the 
orthogonal  projection  of  H2(Z)  in  L2((— joo,  joo);  Z ). 

The  results  depend  on  the  Popov  function  n  defined  by 


n(jcj)  =  G(juj)  +  G(juj)*. 


(2.3) 


Theorem  2.1  Suppose  that  A  is  the  infinitesimal  generator  of  an  exponentially  stable  Co¬ 
semigroup  on  X,  B  E  C{(Cm,X)  and  C  €  C(X,  (Um). 


(a)  IfH(jLo)  >  0  and  n(jo;)  has  invertible  values  for  all  uj  €  M  satisfying 

/OO 

log+  l|n(j<^)  II  i2(  (Z7 7TJ) 


1  +  co2 


-duo  <  oo 


and 


/OO 

iog+  mju) 1  Ht(  (F m) 


du>  <  oo, 

1  +  ur 

then  there  exists  a  spectral  factor  W  E  H'X(C((U  m))  such  that 

n(jw)  =  W(ju)*W(ju)  for  lo  E  R 


(2.4) 

(2.5) 

(2.6) 


(b)  If  there  exists  a  Cw  £  £(D(A),  (Pm)  such  that  for  x  E  X 

Cw(sl  —  A)~1x  E  H2(Wm) 
and  (see  (2.16)  Lemma  2.1  of  Curtain  [5]) 

Ph2  {w{sycw{sl  -  4t)_1x}  =  C{sl  -  A)~1x 

where  VF(s)~  =  W {— s)*;  then 

W{s )  =  Cw(sl-A )“1B, 


(2.7) 


(2.8) 


(2.9) 
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and  for  x  G  D{A) 


A*Qx  +  QAx  =  —C^yCwx 
B*Qx  =  Cx1 


(2.10) 

(2.11) 


where  Q  G  C(X)  is  the  observability  gramian  of  {A,B,Cw)- 

In  general,  (2.7)  and  (2.8)  are  very  difficult  to  check,  but  in  the  case  that  A  is  a  Riesz-spectral 
operator  we  have  the  following  verifiable  conditions. 

Lemma  2.3  Suppose  that  A  has  compact  resolvent  with  eigenvalues  {A n,n  G  IN  }  and 
suppose  that  the  eigenvectors  {en,n  G  IV  }  form  a  Riesz  basis  for  X  and  supn£jy  {i?eAn}  < 
0.  The  following  are  sufficient  conditions  for  the  constrained  Lyapunov  equations  (2.10), 
(2.11)  to  have  solutions  Q  G  C{X)  and  Cw  £  C(D(A)1  X) . 

(*) 

Cl  =  Cwen  =  (1 W(-\n)-l)*Cen]  (2.12) 

(ii)  Either  A  generates  a  holomorphic  semigroup  or  there  exist  numbers  0  <  a  <  b  and 
a  >  0  such  that 

a\Im  Afc|“  <  —Re  A&  <  b\Im  Afc|“  ;  (2-13) 

{in)  there  exists  M  >  0  such  that  for  h  +  iw  G  (T  + 

E  Cw(Cir  <  MR  (2.14) 

-a  neR{h,co)  C,{(Lm) 

where  R{h,uj)  =  {z  G  (T\0  <  Rez  <  h,w  —  h  <  Imz  <  lo  +  h}. 


Conditions  {ii)  and  {in)  ensure  that  C\y  is  an  infinite-time  admissible  observation  oper¬ 
ator  for  T(t).  in  the  terminology  of  “well-posed”  linear  system  (see  Hansen  and  Weiss  [11]). 
Lemma  2.3  covers  a  large  class  of  partial  differential  systems.  We  apply  it  to  our  previous 
Examples  2. A  and  2.B,  both  of  which  have  the  same  structure.  The  transfer  functions  are 


of  the  form 


9{s)  =  E 

n 


(ci  4>n){f>t  4’n) 

S  +  An 


where  bn  =  (5,  <fn)  and  cn  =  (c,  <fn)  are  both  of  order  A  for  large  n  and  Xr. 
examples.  Consequently,  the  analysis  and  the  conclusions  are  the  same. 
Now 


— n2  in  both 


n(jw)  =  E 


+  A  n 


and  since  cn  and  b„  are  of  the  order  of  A  for  large  n.  we  can  estimate  I\-(jio)  ~  -At-  and  it 
is  not  difficult  to  verify  that  (2.4)  holds  and  n  has  a  spectral  factorization  ((2.5)  follows 

from  (2.4)  since  n(jo;)  takes  scalar  values).  In  Curtain  [5]  the  estimate  W{n2)  >  COI);s^ 

nJ 

was  obtained  and  so  |  CR  \  <  const  n~±.  Now  for  our  example,  Cw  will  be  bounded  if 
it  is  of  the  order  n^+e  for  some  positive  e  and  the  admissibility  condition  (2.14)  will  be 
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satisfied  if  \Cyy\  <  const-v/n.  So  G\y  is  unbounded,  but  admissible.  Hence  we  have  satisfied 
all  the  conditions  of  Lemma  2.3  and  there  exist  Q  and  C\y  E  C,(D(A).  (D  m )  satisfying  the 
constrained  Lyapunov  equations. 

We  remark  that,  using  Lemma  2.1  or  Lemma  2.2  on  Example  2. A,  we  obtain  two  dif¬ 
ferent  Qi’s  E  C(X)  satisfying  the  constrained  Lyapunov  equation  (2.10),  (2.11).  Using 
Lemma  2.2  on  Example  2.B  we  also  obtain  a  different  Q 2  E  C(X)  and  L2  E  C(D(A),X) 
satisfying  the  constrained  Lyapunov  equation  (2.10),  (2.11).  In  fact,  for  these  two  examples 
Q 1  and  Q2  are  boundedly  invertible  and  in  section  3  we  shall  see  that  they  are  actually  a 
better  choice  for  our  adaptive  observer  application. 

As  already  noted,  Lemmas  2.1  and  2.2  only  apply  to  special  classes  of  SISO  systems. 
Lemma  2.3  applies  to  a  much  wider  class  of  partial  differential  systems.  Unfortunately, 
Lemma  2.3  is  not  applicable  to  delay  systems.  In  the  following  example  we  show  how  (2.7) 
and  (2.8)  can  be  verified  directly  for  a  delay  system. 

Example  2.C  Consider  the  delay  system 


x(t)  =  —ax(t)  —  bx(t  —  1)  +  u(t)]  a,  b  >  0 
y(t)  =  x(t) 


with  the  transfer  function 


g(s)  = 


1 


s  +  a  +  be 


E  H no. 


Now 


n(jw)  =  g(ju)+g{ju)* 

2  (a  +  6cos(o;)) 


(a  +  6cos(a;))2  +  (lo  —  6sin(u;))2} 
>  0  if  a  >  |6|. 

In  this  case,  it  is  easy  to  find  the  spectral  factor 


W(s)  = 


a  +  f3e 


s  +  a  +  be  s 
W  E  IIoc  and  the  candidate  for  C\,y  is 


a2  +  /32  =  2a,  a/3  =  b. 


(2.15) 

(2.16) 


(2.17) 


{Cw%)  ( t )  =  ax(t)  +  (3x(t,  —  1).  (2-18) 

The  delay  system  (2.15),  (2.16)  can  be  formulated  on  the  state-space  X  =  (D  0  I/2(— 1,0) 
with  generating  operators  defined  by 


Bu  = 


u 


C 


/(•) 


=  r, 


(2.19) 


D(A)  =  { 


/(•) 

df 


E  X\f  is  absolutely  continuous, 


—  (•)  E  L2(-1,0)  and  /( 0)  =  r 


A 


/(•) 


-ar  -  bf(-l) 

d£ 

dO 


(2.20) 
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(see  Curtain  and  Zwart  [8],  chapter  2.4).  We  recall  from  Infante  and  Walker  [12]  that  A 
generates  an  exponentially  stable  semigroup  if  a  —  |5|  >  g  >  0  for  some  positive  constant 
//.  While  C  and  D  are  bounded  operators,  C\y  is  not.  It  is  known  that  it  is  admissible 
(Salamon,  [20]),  which  is  equivalent  to  satisfying  (2.7),  but  we  shall  verify  (2.7)  directly. 
The  resolvent  is  given  by  (Curtain  and  Zwart,  [8],  Lemma  2.4.5) 

)  =  (*?)■  (2'21) 

where 

g(9)  =  esdg( 0)  -  f  es^f(g)dg 
Jo 

d{ °)  =  {r~bJ1  e~s{"+1)f(a) dpj  ■ 

So 

j  =  Cw  ^  j  =  (ag(0)  +Pg(- 1)) 


a  +  (3e  s  (  b(a  +  (3e  s)\ 

A(a)  r+(  A^j- ) 


=  W(s)r  +  (1  -  bW(s))e 


9f{0)d$. 


Now  W  G  and  e“s  G  H0 0  and  so  Cw(sl  -  A)”1  ^  ^  J  G  H2  if  /°x  e~sef(e)d6  G  H2 

and  this  is  a  consequence  of  the  Paley-Wiener  theorem:  J<^1e~Ju'0  f(9)d9  is  the  Fourier 
transform  of  f{9) ![_]., #]• 

Next  we  verify  (2.8).  Since  the  eigenvectors  of  A  span  X  (see  Curtain  and  Zwart  [8] 
Theorem  2.5.10),  it  suffices  to  verify  (2.8)  for  each  eigenvector 


corresponding  to  the  eigenvalue  An: 

A(A„)  =  (An  +  a  +  be~Xn)  =  0  (2.23) 

(see  Curtain  and  Zwart  [8],  Theorem  2.4.6).  Now 

(si  -  A)-1  cf>n  =  — \ —cf>n  and  C(sl  -  A)~X(j)n  = — . 

S  An  S  An 

Next 

cv  4- 

W(srcw(sl  -  A)-1^  =  yCW 

(a  +  /3es)(a  +  /3e~Xn)  1 

A(-s)  s  -  \n 
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Now  the  only  pole  of  W (s)~C\y{sI  —  A)  1</n  in  Res  <  0  is  s  —  Xn  and  so 


PhAW{s)~Cw{sI  ~ 


(a  +  f3eXn )  (a  +  fie  An)  1 

A(— An)  s  —  Xn 

a2  + /32  +  a/3(eXn  +  e~Xn)  1 
— Xn  +  a  +  beXn  s  —  Xn 
2  a  +  b(eXn  +  e~Xn)  1 
2 a  +  b(eXn  +  e~Xn)  s  —  Xn 
from  (2.17)  and  (2.23) 

C(sl  —  A)~x 4>n  as  desired. 


In  fact,  it  is  easily  verified  that  the  solution  to  the  constrained  Lyapunov  equation  (2.10), 
(2.11)  is 

A)- 

For  the  general  retarded  system  with  vector- valued  inputs  and  outputs,  see  Curtain,  [5]. 


3  An  Adaptive  Observer:  main  results 

The  proposed  state  estimator  is 

x(t)  =  A0x(t)  +  Bu(t)  +  BY  (t)y(t)  +  /(f),  2(0)  =  20,  (3.1) 

where  x(t)  is  the  state  estimate  at  time  t  and  T  (t)  is  the  adaptive  estimate  of  the  unknown 
gain.  In  order  to  extract  the  adaptation  rule  for  T(f),  we  use  the  Lyapunov  redesign  method 
[14,  15]  that  has  proved  successful  for  finite-dimensional  systems.  In  this  section,  we  show 
that  the  same  adaptive  observer  that  was  proposed  in  Curtain,  Demetriou  and  Ito  [6]  for 
scalar  (SISO)  systems  can  be  extended  to  the  larger  class  of  multivariable  (MIMO)  systems 
considered  in  this  paper. 

Theorem  3.1  Consider  the  structurally  perturbed  system  (l.l)-(l.S),  where  Ao  is  the  gen¬ 
erator  of  an  exponentially  stable  semigroup  on  X,  B  €  £((Fm,X);  C  E  C(X ,  Cm),  f(t)  is 
a  known  exogenous  signal  which  is  locally  Bochner  integrable,  and  T  is  an  unknown  matrix 
feedback  gain.  If  there  exist  Q  E  C(X)  and  L  E  C(D{A),  X)  or  C(D(A)1  Cm)  satisfying  the 
constrained  Lyapunov  equation  for  x  E  D{A) 

(A0  +  fj.1)*  Qx  +  Q  (A0  +  fj.1)  x  = -L*Lx  (3.2) 

B*Qx  =  CT,  (3.3) 

then  the  state  estimator  defined  by  (3.1)  and  the  adaptation  rule  with  adaptation  matrix 
gain  G  =  GT  >  0  given  by 

f  (f)  =  GCe(t)yT(t) 

f(0)=f0 

have  the  following  properties: 
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(i)  If  u,y  E  1/00(0,00;  (Um),  then  r(t),  xft)  and  the  estimation  error  e(t)  =  x(t)  —  xft) 
are  bounded  in  norm  for  t>  0  and  \\e‘^tQ^  (xft)  —  xft))  \\x  — >  0  as  t  — >  00  ; 

(ii)  If  y  E  L0 o(0,  00;  (C  m)  fl  1/2(0,  00;  W  m)  then  ||e2*  ( x(t )  —  x(t))  ||x  — >  0  as  t  — >  00  ; 


(in)  If  the  conditions  in  (ii)  hold  and,  in  addition,  the  plant  is  persistently  exciting,  i.e., 
there  exist  To,  80  and  eo  such  that  for  each  sufficiently  large  t  >  0  there  exists  t  E 
[ t ,  t  +  Tq]  such  that 


>  eg 


for  every  unit  vector  w  E  Rm ,  then  we  can  achieve  parameter  convergence,  i.e., 


T(t)  -*  r,  as  t  — >■  00. 


Proof  (i)  Consider  the  dynamics  of  the  state  error 

e(t)  =  A0e(t)  +  BTy(t)  -  BT(t)y(t)  ^ 

=  A0e(t)  +  BT(t)y(t) 

e(0)  =  x0-x0  =  eo-  (3-6) 

We  propose  the  following  dynamics  for  the  parameter  error  T(t)  =  F  —  T  (/;) 

f(()  =  —GCe(t)yT{t), 

f  (o)  =  r  -  r.  =  f  0.  M 


First  we  need  to  examine  the  well-posedness  of  the  coupled  system  (3.5),  (3.7)  which  is, 
in  fact,  a  linear  time-dependent  system 


d 

eft) 

dt 

m 

Ao  B[  •  ]y(t) 

-GC(-]yT(t)  0 


e(t) 

m 


The  perturbation  term  D(t, )  :I0(C  rnxm  x  0  C  mxm  is  given  by 

D(t)  = 


0  B[-]y(t) 

~ GC[-]yT(t )  0 


So,  if  y  E  l/oo(0,ti;  (Dm),  (3.8)  has  a  unique  solution  given  by 

=  U(t,  0) 


e(t) 

m 


eo 

f0 


(3.8) 


(3.9) 


(3.10) 


where  U(t,  s )  is  a  mild  evolution  operator  [7]  defined  for  0  <  s  <  t  <  t\.  In  fact,  y(t)  defined 
by  (1.1)-(1.5)  will  always  be  in  (7(0,  ti;(Dm)  for  u,  f  E  Lp(0,  t\;  (Dm),  p  —  1,2  or  00.  In 
general,  eft)  will  not  be  in  D(Aq),  even  if  eo  E  D(Ao).  Sufficient  conditions  for  eft)  E  D(Aq) 
are  that  yft)  E  C1  (0.  t-i-  (Dm),  which  are  very  strong.  However,  we  assume  this  initially  to 
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facilitate  the  Lyapunov  argument.  We  examine  the  asymptotic  properties  of  (3.10)  using 
the  following  Lyapunov  functional  for  (  ~  ^ 


V(e,T )  =  (e,Qe)+Tr{rTG"1r},  (3.11) 

where  Q  is  the  solution  to  (3.2).  Since  e(t)  E  D(Aq)  we  may  differentiate  V  along  solutions 
of  (3.8)  for  0  <  t  <  t\  to  obtain 

V (e,  f )  =  (A^e  +  BTy,  Qe)  +  (Qe,  A$e  +  BTy)  +  2Tr  |fTG_1f  | 

=  —  ||Le||2  —  2fj(e,Qe)  —  2{Ce)TYy  +  2Tr  {y(C'e)Tf }  (3.12) 

using  Lemma  2.1  (3.2),  (3.3)  and  (3.7) 

=  —  ||Le||2  —  2y(e,Qe)  using  (2.9)  and  bT a  =  Tr(abT). 


We  now  integrate  (3.12)  from  t  =  0  to  t  =  t\  to  obtain 

rh 


i(ti),Qe(ti))  +Trjrr(fi)G  1r(<i)|  +  ^  \\Le(t)\\2  dt  +  2/j  (Qe(t) , e(t))  dt 

=  (eo,Qe0)+Tr{f^G-1fo}. 


(3.13) 


Notice  that  although  we  have  assumed  that  eo  E  D(Aq)  and  y  E  CL(().  /y:  d'm)  to  derive 
(3.13),  all  terms  make  perfectly  good  sense  for  eo  E  X  and  y  E  (7(0, fi;  € m ) .  Moreover, 
(3.10)  and  the  facts  that  sup0<s<t<t1  \\U(t,  s)||  <  oo,  and  that  D{Aq)  is  dense  in  X,  show 
that  (3.13)  can  be  extended  to  all  eo  E  X.  We  now  extend  (3.13)  to  all  y  E  Loo (0,2;  ®m) 
by  appealing  to  Lemma  A.l  in  the  appendix,  which  shows  that  if  we  approximate  y  by  a 
sequence  yn  E  C71  (0,  ;  <D  ”*)  satisfying 


[  \\y(s) -yn{s)\\2ds 

Jo 


0,  as  n  — >■  oo, 


then  there  holds 


sup  \\U(t,  s)  —  Un(t,  s)  ||  — >  0  as  n  — >  oo. 

0<s<t<ti 


(3.14) 


So  the  respective  solutions  to  (3.8)  satisfy 


sup  (||e(2)  -  en(t)\\x  +  ||r(2)  -  Yr 

0<s<t<ti  v 


—>•0  as  n  oo 


(3.15) 


and  this  suffices  to  show  that  (3.13)  holds  for  any  y  E  Loo(0, 2i;  (Dm)  and  eo  E  X.  This 
implies  that  f  €  Loo(0,  oo;  <D  mxm ). 

Next,  we  define  q(t)  :=  ||Q2e(t)||2  and  deduce  the  following  from  (3.13) 

q{t i)  +  2 n  £  f{s)ds  <  /( 0)  +  Tr  {r^G'-1r0}  .  (3.16) 

Equation  (3.16)  implies  that  q(ti)  <  e_2/itlE(0)  or  equivalently 

||e^Q^e(2i)||2<  HQ^eoll2.  (3.17) 
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Now  t\  can  be  chosen  arbitrarily  large  and  so  e(t)\\2  — >  0  as  f  — >  oo. 

Finally,  since  T(t)  and  y(t)  are  uniformly  bounded  in  norm  for  t  >  0  and  Aq  generates 
an  exponentially  stable  semigroup,  (3.5)  shows  that  e(t)  is  uniformly  bounded  in  norm  for 
t  >  0.  Similarly,  (2.3)  shows  that  x(t)  is  uniformly  bounded  in  norm  for  u  E  Foo(0,  oo;  Cm) 
and  hence  x(t)  =  e(t)  +  x(t)  has  the  same  property. 

( ii )  If  y  E  1/2(0,  oo;  (D m),  then  from  (i),  the  forcing  term  BT(t)y(t)  in  (3.5)  is  in 
L2(0,  oo]  X).  This,  together  with  the  fact  that  Ao  generates  an  exponentially  stable  semi¬ 
group  implies  that  e(t)  — »  0  as  t  — »  oo  (see  Lemma  12  in  Curtain  and  Oostveen  [16]). 

(in)  The  parameter  convergence  is  proven  by  applying  the  results  in  Baumeister  et  al  in 
[2]  to  equation  (3.8).  We  discuss  briefly  the  minor  subtleties  required  to  apply  Definitions  2.1 
and  3.3  and  Theorem  3.4  of  [2],  The  definition  of  the  (admissible)  plant  requires  that 
\(BTy(t),ip)x\  <  q(y) ||r|| \ip\x-  Since  the  operator  B  is  bounded  and  the  output  y  E 
//^(O,  oo;  (D  m)  we  have  that  indeed  the  plant  satisfies  Definition  2.1  of  [2]  with  q(y)  = 
l\y\Loo-  The  persistence  of  excitation  condition  in  Definition  3.3  of  [2]  takes  the  form  of 


sup 

IfolKl 


/•f+<5o 

J  (BTy(T),ip)di 


>  eo, 


for  all  T  EHmxm  with  ||r||  =  1,  which  simplifies  to 


sup 

WIKI 


rt+6  o 

J_  yT(r)rT(B»dr 


This  finally  becomes 


rt+So 

j  y  (r)wdr 


>  e0- 


for  some  w  €  IRm  with  unit  norm.  Then  a  direct  application  of  Theorem  3.4  in  [2]  yields 
the  required  result.  This  completes  the  proof  of  Theorem  3.1.  □ 


4  Adaptive  Compensators 

In  this  section,  we  propose  an  adaptive  compensator  for  the  perturbed  plant  (1.1)-(1.3) 
where  f(t)  is  a  known  exogenous  signal.  We  obtain  results  for  f(t)  a  periodic  signal  and 
for  /  E  1/2(0,  oo;  X).  First  we  apply  output  injection  to  obtain  a  modified  control  problem: 

u(t)  =  u2-V(t)y(t).  (4.1) 

This  has  the  advantage  of  producing  the  new  estimator  dynamics 

x(t)  =  A0x(t)  +Bu2(t)  +f(t) 
a;(0)  =x0 


and  the  same  error  dynamics  (3.8)  for 


e(t) 

m 


as  before. 


So  it  remains  to  design  a  controller  U2(t)  for  the  system  (4.2). 
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We  use  the  LQR  control  design  from  Da  Prato  and  Ichikawa  [19].  Suppose  that  (To,  B,  C2) 
is  exponentially  stabilizable  and  detectable  and  R  =  R1  >  0  G  £(C’ m).  We  seek  to  minimize 
the  average  cost 

J^  =  T^ol'Io  (\\C2x(t)\\2 +  \\R-^u{t)\\2^dt  (4.3) 

over  all  controls  satisfying  limj^oo  ^  Jq  \\u(t)\\2dt  <  00  and  for  which  the  corresponding 
closed  loop  trajectory  is  bounded  on  t  >  0.  They  showed  that  if  f(t )  is  periodic,  the 
minimizing  control  law  is  given  by 

u2(t)  =  —R~1B*  ( Px(t )  +  r(i)) ,  (4.4) 

where  P  =  P*  G  C(X)  is  the  solution  to  the  Riccati  equation  for  x  G  D(Aq) 

AqPx  +  PA0x  -  PBR~1B*Px  +  C2C2x  =  0  (4.5) 

and  r(t)  is  the  solution  to 

r(t)  =  —  PBR~1B *)  r(t)  —  Pf  (t),  r(t)  — >  0  as  t  — >■  00.  (4.6) 

Equation  (4.6)  has  the  solution 

/OO 

Tp{s  —  t)P  f  (s)ds,  (4.7) 

where  Tp(t)  is  the  exponentially  stable  Co-semigroup  generated  by  A  —  BR~1B*P. 

The  closed  loop  trajectory  converges  exponentially  fast  to  the  periodic  solution 

p(t)  =  J*  TP(t- s)  (/(a)  -  BR-'B'risj)  ds ,  (4.8) 

i.e., 

lim  evt\\x(tp0)  -p(t)\\x  =  0,  (4.9) 

r— >■  00 

where  u  is  the  decay  rate  of  Tp(t). 

In  the  case  of  a  constant  exogenous  signal,  f(t)  =  /o, 

p(t)  =  -(au-  BR~1B*Py1A*0  (Rq  -  BR~1B*P)  f0. 

We  note  that  if  /  G  1/2(0,  00;  X),  then  the  feedback  control  law 

u2(t)  =  -il_1B*Pa(<) 

has  the  effect  of  ensuring  that  e^x^t)  — >■  0  as  t  — >  00  (see  Lemma  12  in  Curtain  and 
Oostveen  [16]). 

We  propose  the  following  adaptive  compensator  for  the  case  of  a  known  periodic  exoge¬ 
nous  input 

x(t)  =  ^4o  —  BR~1B*P^J  x(t)  —  BR~1B*r(t )  +  /(f);  x(0)  =  xq 

r(t)  =  (To  -  PRiT1#*)  r(t)  -  Pf{t),  r(t )  ->  0  as  t  ->  00 
u{t)  =  —R~1B*  ( Px{t )  +  r(t))  -  f  (<)i/(<) 
f  (<)  =  GCe(t)yT(t),  f  (0)  =  P0 
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for  our  structurally  perturbed  plant  (1.1)-(1.3). 

In  section  3  we  showed  that 

||e24Q2  (x(t)  —  x(t))  \\x  — >  0  as  t  — »  oo 

independently  of  the  choice  of  the  control.  So  combining  this  with  the  results  in  this  section, 
we  conclude  that  for  the  case  of  a  known  periodic  input  f(t) 

|| (■ x(t )  —  x(t))  ||x  — >  0  as  t  -»  oo, 

where  (3  =  min(zA  and  p(t)  is  given  by  (4.8).  If  f(t)  is  in  L2(0,  oo;X),  then 

We&Q*  (x(t)  —  x(t))  \\x  ->  0  as  t  — »  oo. 


5  Examples  and  Numerical  Results 


We  present  some  numerical  results  for  the  three  examples  considered  in  section  2.  For  each 
of  these  examples,  there  exists  a  solution  Q  6  C(X)  satisfying  (3.2)  for  a  certain  p  >  0, 
and  in  all  three  cases  Q  is  invertible.  Consequently,  we  can  conclude  that  for  the  adaptive 
observer  (3.1)  with  adaptation  rule  (3.4) 

||ea*  ( x(t )  —  x(t))  llx  — >  0 

and  with  the  adaptive  compensator  of  section  4 

\\em  (x(t)  - p(t ))  \\x  ->  0. 


All  the  computations  described  below  were  carried  out  on  a  Digital  Personal  Work¬ 
station  433  a u-Series  in  the  Mechanical  Engineering  Department  at  Worcester  Polytechnic 
Institute.  A  finite  element  Galerkin  approximation  scheme  based  on  spline  elements  was 
used  for  the  spatial  discretization  of  the  two  PDE’s  similar  to  the  one  developed  in  [2], 
The  resulting  finite  dimensional  ODE  systems  were  integrated  in  time  using  a  Fehlberg 
fourth-fifth  Runge-Kutta  method.  The  delay  system  in  Example  2.C  was  discretized  using 
the  method  presented  in  the  paper  by  Ito  and  Kappel,  [13].  The  resulting  evolution  (finite 
dimensional)  system  was  similarly  integrated  using  the  Runge-Kutta  code  rkf45.f. 
Example  2. A  As  was  already  mentioned  in  section  2,  we  can  choose  in  this  case  p  =  ir2  —  e 
in  (2.1)  and  define  the  operators  Q  and  L  via  (2.2).  Alternatively,  when  Lemma  2.1  is  used, 
we  have  Q  =  I  and  L  =  0  with  the  same  p.  The  input  operator  b(x)  was  chosen  as 


b(x) 


1  on  [0, 1) 

0  elsewhere 


The  unknown  gain  was  chose  as  T  =  1  and  as  initial  conditions  we  chose  z(x,  0)  =  sin(7r.x) 
and  z(x,  0)  =  cos(2-7ra:)  —  1.  The  exogenous  input  f(x,t)  =  50x[o,i]  {x)  sin(27rf).  The  initial 
guess  for  T(0)  =0  with  an  adaptive  gain  of  G  =  20.  Figure  1(a)  depicts  the  time  evolution 
of  the  output  state  error  Ce(t)  =  Jq(z(x.  i)  —z(x.  t))dx.  The  convergence  to  zero  is  achieved 
within  0.5  seconds.  The  parameter  estimate  T(f)  (dashed)  and  the  actual  value  of  V  =  1 
are  depicted  in  Figure  1(b).  Parameter  convergence  is  achieved  in  4  seconds. 
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Evolution  of  output  error,  Ce(t) 


Time  (sec) 


Figure  1:  Evolution  of  (a)  output  error  and  (b)  parameter  estimate  T(f)  (dashed)  -  actual 
parameter  T  (solid). 


Example  2.B  Equations  (2.1)  and  (2.2)  can  be  satisfied  with  //  =  ^-n2—e  where  the  param¬ 
eter  a  =  0.2.  Initial  conditions  were  set  as  z(x,  0)  =  sin(-7nr)  and  z(x,  0)  =  —  0.25  sin(7nr).  A 
constant  in  space  and  time  exogenous  function  is  implemented  as  f(x,t)  =  50x[o,i](a:)  and 
T(0)  =  0  with  G  =  2.  With  these  values  of  initial  conditions,  it  is  observed  in  Figure  2(a) 
that  the  the  output  state  error  converges  to  zero  in  3.5  seconds.  Furthermore,  the  parameter 
T(f)  converges  to  the  actual  value  T  =  1  in  4  seconds,  as  shown  is  Figure  2(b). 

Example  2.C  The  plant  parameters  were  chosen  as  a  =  3,  b  =  1.  In  this  case  the  solution 


to  the  constrained  Lyapunov  equations  (2.10),  (2.11)  is  Q  = 

The  actual  value  of  the  parameter  was  T  =  0.4  with  initial  condition  for  its  estimate 
chosen  as  T(0)  =  0.2.  The  initial  state  was  set  at  x(t  —  1)  =  sin(4f  —  1)  —  sin(— 1)  and  the 
state  estimate  as  x(t  —  1)  =  0.5sin(4f  —  1);  thus  x(0)  =  sin(3)  +  sin(l)  and  x(0)  =  0.5sin(3). 
Here  we  had  /(f)  =  0  for  the  exogenous  signal  and  chose  an  adaptive  gain  of  G  =  500.  It  is 
observed  from  Figure  3(a)  that  the  state  estimate  converges  to  the  plant  state  in  about  2 
seconds.  Parameter  convergence  is  also  achieved  at  about  5  seconds.  For  numerical  results 
for  a  multivariable  example  the  reader  is  directed  to  [9] . 


I  0 
0  (3  ±  y/8 )I 


References 

[1]  A.  V.  Balakrishnan,  On  a  generalization  of  the  Kalman- Yakubovic  lemma ,  Applied 
Mathematics  and  Optimization,  31  (1995),  pp.  177-187. 


16 


Evolution  of  output  error,  Ce(t) 


Time  (sec) 

Figure  2:  Evolution  of  (a)  output  error  and  (b)  parameter  estimate  r (/;)  (dashed)  -  actual 
parameter  T  (solid). 


[2]  J.  Baumeister,  W.  Scondo,  M.  Demetriou,  and  I.  Rosen,  On-line  parameter 
estimation  for  infinite  dimensional  dynamical  systems ,  SIAM  J.  Control  and  Optimiza¬ 
tion,  3  (1997),  pp.  678-713. 

[3]  R.  F.  Curtain,  Corrections  to  The  Kalman- Yakubovich- Popov  Lemma  for  Pritchard- 
Salamon  Systems ,  Systems  and  Control  Letters,  28  (1996),  pp.  237-238. 

[4]  - ,  The  Kalman-  Yakubovich- Popov  Lemma  for  Pritchard- Salomon  Systems ,  Systems 

and  Control  Letters,  27  (1996),  pp.  67-72. 

[5]  - ,  Linear  operator  inequalities  for  stable  weakly  regular  linear  systems ,  submitted, 

(1997). 

[6]  R.  F.  Curtain,  M.  A.  Demetriou,  and  K.  Ito,  Adaptive  observers  for  structurally 
perturbed  infinite  dimensional  systems ,  in  the  Proceedings  of  the  1997  IEEE  Conference 
on  Decision  and  Control,  San  Diego,  California,  1997. 

[7]  R.  F.  Curtain  and  A.  J.  Pritchard,  Infinite  Dimensional  Linear  Systems  Theory , 
Lecture  Notes  in  Control  and  Information  Sciences,  Vol.  8,  Springer- Verlag,  Berlin, 
1978. 

[8]  R.  F.  Curtain  and  H.  J.  Zwart,  An  Introduction  to  Infinite  Dimensional  Linear 
Systems  Theory ,  Texts  in  Applied  Mathematics,  Vol.  21,  Springer- Verlag,  Berlin,  1995. 


17 


Evolution  of  x(t)  and  x(t) 


Figure  3:  Evolution  of  (a)  plant  state  x(t)  (solid)  and  state  estimate  x(t)  (dashed);  (b) 

parameter  estimate  T(f)  (dashed)  -  actual  parameter  T  (solid). 

[9]  M.  A.  Demetriou,  R.  F.  Curtain,  and  K.  Ito,  Adaptive  observers  for  structurally 
perturbed  positive  real  delay  systems ,  in  the  Proceedings  of  The  Fourth  International 
Conference  on  Optimization:  Techniques  and  Applications,  Curtin  University  of  Tech¬ 
nology,  Perth,  Australia,  July  1-3  1998. 

[10]  M.  A.  Demetriou  and  K.  Ito,  Adaptive  observers  for  a  class  of  infinite  dimensional 
systems ,  in  Proceedings  of  the  13th  World  Congress,  IFAC,  San  Francisco,  CA,  July 
1996. 

[11]  S.  Hansen  and  G.  Weiss,  New  results  on  the  operator  Carleson  measure  criterion , 
IMA  Journal  of  Math.  Control  and  Information,  14  (1997),  pp.  3-32. 

[12]  E.  Infante  and  J.  Walker,  A  Lyapunov  functional  for  scalar  differential  difference 
equations ,  Proc.  Royal  Society  of  Edinburgh,  79A  (1977),  pp.  307-316. 

[13]  K.  Ito  and  F.  K appel,  A  uniformly  differentiable  approximation  scheme  for  delay 
systems  using  splines ,  Applied  Mathematics  and  Optimization,  23  (1991),  pp.  217-262. 

[14]  H.  K.  Khalil,  Nonlinear  Systems ,  Macmillan,  New  York,  1992. 

[15]  K.  S.  Narendra  and  A.  M.  Annaswamy,  Stable  Adaptive  Systems ,  Prentice  Hall, 
Englewood  Cliffs,  NJ,  1989. 

[16]  J.  C.  Oostveen  and  R.  F.  Curtain,  Riccati  equations  for  strongly  stabilizable 
bounded  linear  systems ,  Automatica,  to  appear,  (1998). 


18 


[17]  L.  Pandolfi,  The  Kalman-  Yacubovich  theorem:  An  overview  and  new  results  for 
hyperbolic  boundary  control  systems ,  Rapport  Interno  15,  Dipartimento  di  Matematica, 
Politecnico  di  Torino,  Torino,  1996. 

[18]  - ,  Dissipativity  and  Lur’e  problem  for  parabolic  boundary  control  systems ,  Rapport 

Interno  1,  Dipartimento  di  Matematica,  Politecnico  di  Torino,  Torino,  1997. 

[19]  G.  D.  Prato  and  A.  Ichikawa,  Quadratic  control  for  linear  periodic  systems ,  Appl. 
Math.  Optim.,  18  (1988),  pp.  39-66. 

[20]  D.  Salamon,  Control  and  observation  of  neutral  systems ,  in  Research  Notes  in  Math¬ 
ematics,  91,  Pitman  Advanced  Publishing  Program,  Boston,  1984. 

[21]  O.  J.  Staffans,  Quadratic  optimal  control  of  stable  abstract  linear  systems ,  in  IFIP 
conference,  Modelling  and  Optimization  of  DPS  with  Applications  to  Engineering, 
Warsaw,  Poland,  1995. 

[22]  - ,  Coprime  factorizations  for  well-posed  linear  systems ,  SIAM  Journal  of  Control 

and  Optimization,  (1997),  p.  to  appear. 

[23]  - ,  Quadratic  optimal  control  of  stable  well-posed  linear  systems ,  Transactions  of  the 

American  Mathematical  Society,  (1997),  p.  to  appear. 

[24]  - ,  Quadratic  optimal  control  of  well-posed  linear  systems ,  SIAM  Journal  of  Control 

and  Optimization,  (1997),  p.  to  appear. 

[25]  M.  Weiss,  Riccati  Equations  in  Hilbert  Spaces:  A  Popov  function  approach ,  PhD 
thesis,  Rijksuniversiteit  Groningen,  The  Netherlands,  1994. 

[26]  - ,  Riccati  equation  theory  for  Pritchard-Salamon  systems:  a  Popov  function  ap¬ 

proach,  IMA  Journal  of  Math.  Control  and  Information,  14  (1997),  pp.  45-83. 

[27]  M.  Weiss  and  G.  Weiss,  Optimal  control  of  stable  weakly  regular  linear  systems, 
Mathematics  of  Control,  Signals  and  Systems,  10  (1997),  pp.  287-330. 


19 


Appendix 


Lemma  A.l  Suppose  that  A  generates  a  Co-semigroup  on  the  Hilbert  space  X  and  consider 
the  mild  evolution  operator  U(t,s)  generated  by  A  +  Ei=i  A;  £  L(X)  and  tji  € 

Loo(0,<i).  Let  Un(t,s)  be  the  evolution  operator  generated  by  A  +  where  for 

each  i  yf{t)  is  a  sequence  of  functions  in  C1( 0,ti)  satisfying 

rt  i 

/  WViit)  ~  yi(t)\\2dt  -»  0  os  n  -M». 

Jo 

There  holds 

sup  ||[/(t,  s)  —  [/„(<,  s)||£(x)  — t  °°  asn-^oo. 

0<s<t<t\ 


Proof  We  only  give  a  detailed  proof  for  k  —  1,  since  the  arguments  extend  readily  to 
any  finite  k.  We  recall  from  Curtain  and  Zwart  [8]  the  defining  equations  for  U(t,s)  and 
Un{t,s):  ^ 

U(t,  s)x  =  T(t  —  s)  +  (  T(t  —  a)Dy(a)U(a,s)xda  (A.l) 

J  S 

Un(t,  s)x  —  T(t  —  s)  +  f  T(t-  a)Dyn(a)Un(a,s)xda  (A.2) 

J  S 


and  the  estimate 

||C/(a,  s)  ||  < 

where  \\T(t)  ||  <  Mewt,  t>  0, 


(A.3) 

(A.4) 


and  n  =  M||P||||i/||Loo(0)tl)  >  0. 


Consider  the  following  estimates  obtained  using  (A.l)  and  (A.2) 


\\U(t,s)-Un(t,s)\\ 


<  [  \\T(t  -  a)D\\\y(a)  -  yn(a)\\\U(a,s)\\da 

J  S 

+  [  ||P(*  -  a)D\\\\yn\\Loo\\U(a,s)  -  Un(a,s)\\da 

J  S 

<  ||A||  f  Me^-^Me^+^-^lyia)  -  yn{a)\da 

J  S 

+  11^1111^11^  f  Me“^\\U{a,s)  -  Un(a,  s)\\da. 

J  S 


20 


Defining  fn(t,s)  =  e  s^\\U(t,  s)  —  Un(t,s)\\  we  obtain 

fn(t:s)  <  M2\\D\\  [  eftae~^s\y(a)  -  yn(a)\da 

J  S 

+  11-^1111^11^  /  fn{a,s)da 

J  S 

t  I  t  - 

<  M2\\D\\(^j^  \y{a)  —  yn{a)\2  do^j  Q [  e2^a~s) da'j 

+  PIIIIi/t,IIl00  [  fn(a,s)da 

J  S 

=  C\\e2^  -  l\h  (£  | y(a)  -  yn(a)\2da )  * 

+C2  [  fn(a,s)da 

J  S 

where  C'i  and  C2  only  depend  on  t\. 

Thus 

fn{t,s)  <  2C'ie/i(t_s)||y-yn||i2(0,t1)  +  C2  [  fn{a,s)ds 

J  S 

and  differentiating  with  respect  to  t  for  fixed  s  yields 

^r(M)  <  2Cf/^(t“+||y-yJi2(0;tl)  +  C2fn(t,s) 

and 

f  (e-Ca'/»(t,  »))  <  -  yn ||M„,tl|. 

We  integrate  from  t  to  s  noting  that  fn{s.  s)  —  0  to  obtain 

e~C2tfn{t,s)  <  2C1ye~IJ's  f  e^~C2^d/3\\y  -  yn\\L2^tl) 

J  S 

=  (e<^Cl,‘ “  e<^Cl)')  »»  -  »»IImo,<.) 


and 


/»(«, »)  <  (e'z-c-O*  _  ||„  _  „„||fa(0i(l) 


and 

||C/(t,s)-C/„(t,S)||  < 
which  proves  our  claim. 


2C\y 

n-c2 


g (w+/j>)(t-s)  _  e(w+C2)(t-s ) 


III/  “  2/ralU2(0,ti) 
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